Electron-phonon interaction on bundled structures: static and transport properties 



Ilaria Meccoli 

Dipartimento di Fisica, Universitd di Parma, and Istituto Nazionale per la Fisica della Materia (INFM) 

Parco Area delle Scienze 7a, 1-43100 Parma, Italy 



Unitd di Parma, 



Massimo Caponc 

International School for Advanced Studies (SISSA), and Istituto Nazionale per la Fisica della Materia (INFM) 
Unitd Trieste-SISSA, Via Beirut 2-4, 1-34014 Trieste, Italy 
(February 1, 2008) 



X 



We study the small-polaron problem of a single electron interacting with the lattice for the Holstein 
model in the adiabatic limit on a comb lattice, when the electron-phonon interaction acts only on the 
base sites. The ground state properties can be easily deduced from the ones of a linear chain with an 
appropriate rescaling of the coupling constant. On the other hand, the dynamical properties, that 
involve the complete spectrum of the system, present an "exotic" behavior. In the weak coupling 
limit the Drude weight (zero-frequency conductivity) is enhanced with respect to its free-case value, 
contrary to the linear chain case, where for every finite value one has a suppression of the Drude 
peak. More interestingly, the loss of coherent electron motion and the polaronic localization of the 
carrier occurs for different coupling values. Thus for intermediate coupling, a novel phase appears 
with large kinetic energy and no coherent motion. 



71.23.-k, 71.55.Jv, 71.38.+i, 63.20.Kr 
I. INTRODUCTION 

In the last years, a general interest has been growing 
around the study of physical properties of inhomogeneous 
discrete structures. The general purpose is to character- 
ize the behavior of systems like amorphous solids, glasses, 
polymers and biological systems in general, where one 
can expect the specific geometry to dramatically influ- 
ence physical properties, ranging from transport and dif- 
fusive properties, to the thermal or electrical ones. A 
very general and powerful formalism has been developed, 
the so-called random walk problem, where the structure 
is explored by a classical walker, that randomly jumps 
to an arbitrary nearest-neighbor site at every time step. 
This approach can be thought as the discrete version of 
the classical Boltzmann equation, and the asymptotic be- 
havior of the average distance of the walker from the 
starting site gives the diffusion law for the considered 
structure. On translationally invariant graphs, the diffu- 
sion law predicts that this distance scales with the time 
step with a power law depending on the euclidean dimen- 
sion of the lattice. On more general graphs, the diffusion 
law is governed by a new parameter characterizing the 
structure, generally known as spectral dimension j[p|§, 
that can be different from the euclidean dimension in 
which the graph is eventually embedded, giving rise to 
the so-called anomalous diffusion . One of the main 
difficulties of this approach is that the explicit calcula- 
tion of this important parameter may not be trivial for 
general structures. 

Nevertheless, there exists a wide class of graphs, called 
bundled structures, where explicit calculations can be 
performed. The random walk and the oscillation prob- 



lems have already been studied on these structures 
and an anomalous diffusion law for a classical walker 
moving along the base has been discovered. Further- 
more, the nearest-neighbor tight binding model has been 
explicitly solved for electrons moving on these graphs, by 
means of an exact resummation of the perturbative ex- 
pansion in the hopping parameter Q . This perturbative 
approach becomes very complicated, and in a last in- 
stance useless, when the electrons experience any kind of 
interaction, both between them, and with external fields, 
like e.g. phonons. This work is devoted to the first ex- 
tension to the interacting case of the problem of electrons 
moving on bundled structures. We restrict ourselves to 
the case of a single electron which interacts with local 
oscillators only on the base sites of a comb lattice. This 
is a first step towards the detection of some quantum 
counterpart of the anomalous diffusion laws. 

The subject of electron-phonon (e-ph) interaction on 
unconventional structures has been already addressed 
mainly within the framework of the Su-Schrieffer-Heeger 
H model, originally introduced to describe the properties 
of poliacetylene chains j[o| . We will attack the problem 
of e-ph interaction on unconventional structures from an- 
other point of view. Rather than a realistic description 
of actual compounds, we aim to provide a general picture 
of the mutual effect of e-ph coupling and geometric com- 
plexity. We will therefore focus on the Holstein model 
JrT| , which despite its simplicity, is expected to capture 
all the main physical properties of more general inter- 
actions. In section || we sketch the basic formalism for 
bundled structures, and show how the problem can be 
mapped onto a base-only one. In section III the small- 



polaron problem on the comb lattice is discussed, with 
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particular emphasis on the dynamical properties. Section 
IV is devoted to conclusions and future perspectives. 



II. GENERAL FORMALISM 

A bundled structure can be built joining each point of 
a base graph B with a copy of a fibre graph T in such a 
way that there is only one common site between B and 
T . As previously explained, we shall limit our study to 
the simplified case in which the particles experience in- 
teraction only on the base sites. Within this hypothesis, 
the generic hamiltonian for a bundled structure assumes 
a block form: 



Hb Hbf 
Hfb Hp 



(1) 



Hb and Hp act on the base and the fibres sites respec- 
tively, while Hbf — Hj. B are rectangular matrices con- 
taining the hopping terms joining the fibres to the base. 
Consequently, a generic state can be written as 



I''/) 



(2) 



with belonging to the base and \rf) to the fibres. Since 
this separation of degrees of freedom is completely for- 
mal, we must ensure that it does not introduce unphysical 
states in the Hilbert space. So we impose a normaliza- 
tion constraint: (ip\ip) + (i]\r]) — (4>\4>) = 1. Let us write 
the eigenvalues equation for this structure: 



H B \ip e ) + H BF \r) e ) = e\ip e ) 
H F \r, € )+H FB \^ e ) =e|^) 



(3) 



Because Hp does not contain any interaction term, one 
can get rid of \ij e } from (||), reducing to a "base-only" 
problem: 



(e-H B - H BF {e - Hp^Hps)^ 
(Ml - H BF(e - H f )- 2 Hfb}\A 







(4) 



Explicitly using the simple form of the rectangular matrix 
Hbf (reflecting the peculiar geometry of the bundled 
structure), it can be shown that the composite operators 
in (Q) have only diagonal terms simply expressed in term 
of the fibre-only Green's functions: 

[H BF (e - Hp)- x Hf B ].. = 6 l3 n f t 2 F 00 {e) (5) 
[H BF {e- Hp)- 2 H FB ].. = 5 l3 n f t 2 Y,F 0m (e)F m0 (e) (6) 

Here «/ is the number of independent fibres joined to a 
single site of the base; t is the hopping constant relative 
to the link joining the fibre and the base; i*; m (e) = [e — 
H^llm * s ^ ne fit> re Green's function , and i — indicates 



the first site of the fibre when coming from the base. So 
we have reduced ourselves to an effective problem for the 
only base sites: 



(e-H B -n f t 2 F 00 (e)) |^ 6 ) = 



l + n f t 2 £ m F 0m (e)F m o(e) 
We can rewrite Eq.(^) as 

(/(e) - H B ) \i> e ) = 0, 



n(e). 



(7) 

(8) 

(9) 



defining /(e) = e — nft 2 F 00 (e). Two effects keep track 
of the presence of the fibres. A time-dependent potential 
V(t) = n f t 2 J di J ip*(T)F a(T - T , )ip(i J ), of purely geo- 
metrical origin, acts on the base, accounting for the pos- 
sibility for the particle to explore the fibres when moving 
on the base. On the other hand, the normalization of the 
eigenvectors n(e) depends self-consistently on the corre- 
sponding eigenvalue. This fact will have a dramatic con- 
sequence in presence of interaction, and the Schrodinger 
equation contains a nonlinear term. Both these aspects 
will be analyzed in the subsequent section. 

Let us note that the decomposition (|], ||) holds under 
the very general hypothesis that the fibres and the base 
are connected only by an hopping term. If Hp contains 
instead an interaction term, we generally not know the 
explicit form of the Green's functions Fij (e) and conse- 
quently the expressions (||), (||) can be not be computed 
exactly. 



III. ELECTRON-PHONON INTERACTION: THE 
COMB LATTICE CASE 

In this section we turn to the problem of electron- 
phonon interaction on bundled structures. We will con- 
sider the simplest model for e-ph interaction, namely the 
Holstein molecular crystal model Jll]]. 



g y^rti(aj + a\) + cj a\ai. 



(10) 



In this model tight-binding electrons interact with local 
dispersionless oscillators of frequency wo, g is an e-ph 
coupling between the displacement a% + a\ of the oscilla- 
tor and the local electronic density n, = J2„ c la c icr- O n 
really general grounds, the model is expected to display 
a polaronic behavior for large el-ph coupling. A polaron 
is a bound state of electrons and phonons in which the 
electron moves carrying with itself a phonon cloud that 
strongly diminishes its mobility, eventually giving rise to 
localization 11 liq]. Since, at least for g = 0, the elec- 
tron motion is completely free, some kind of transition or 
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crossover from a weak-coupling free-carrier state, and a 
strong-coupling polaronic state must occur at some cou- 
pling. Despite the simplicity of the model, even the case 
of a single particle interacting with the set of oscillators 
is a non-trivial many-body problem that can be solved 
only in particular limits. Wc mention the solution ob- 
tained in Ref. [|l3| within the Dynamical Mean Field, 
which is exact in the limit of infinite spatial dimensions, 
and the theorem by Gerlach and Lowen that claims that 
no phase transition occurs for finite luq In general 

cases one has to resort to approximate techniques or nu- 
merical studies p5|-p7| . 

The nature of the transition between a free carrier (for 
g = 0) and a polaron depends on the value of Uo/t and 
spatial dimensionality jl6|. In fact, even for a single par- 
ticle we have two independent parameters in the Hamil- 
tonian (|Io|), which determine different physical regimes. 
In particular, while g is a measure of the strength of the 
e-ph interaction, 7 = u>o/t is a measure of the adiabatic- 
ity of the system, i.e. it controls the relative value of the 
typical phononic and electronic energy scales. It has been 
shown that in the adiabatic regime 7 <C 1, the condition 
that rules the polaronic crossover is A = g 2 /(2ojotd) > 1, 
where d is the dimensionality. This condition is simply 
understood as an energy convenience condition as soon 
as we notice that the polaronic energy is E po i = —g 2 /u)Q, 
and the free carrier energy is Ef ree = —2td. Notice that 
Ef ree contains the total kinetic energy p5|Jl^ |. 

We will restrict ourselves to the extreme adiabatic limit 
u>o/t = 0, in which the lattice degrees of freedom become 
classical. In this limit, the Hamiltonian can be recast in 
the form 



H = 



-*£(< 

(i,i)><7 



ia^]° 



h.c) + g ^ n i x i + Tj^y^/. 



(11) 



Here g = y / 2k/ujog, k is the elastic constant, and Xi are 
the classical variables associated to the local displace- 
ments of the oscillators. It can be shown that the relevant 
adimensional coupling the adiabatic limit is A = g 2 jAkt. 

Kabanov and Mashtakov pq| have solved the model for 
a single particle for 1-2-3 dimensional systems showing 
that in the one-dimensional case the system has a local- 
ized bound state for arbitrarily small e-ph coupling, ev- 
idencing only a crossover between a weak-coupling large 
polaron extended on many lattice sites, and a strong- 
coupling small polaron localized on a single site. In two 
and three spatial dimensions a first-order transition (level 
crossing) occurs between a delocalized state and a small 
polaron one. We will generalize the approach of refer- 
ence [|l8| to the case of a comb lattice, a specific bundled 
structure, in which the base is a linear chain, and two 
semi-infinite chains are linked to each site of the chain 
(fibres). We will also consider only the case in which the 
e-ph interaction is limited to the base. This lattice can 



be viewed as a good model for describing the geometry 
of many branched polymers ]l9| , and, on the other side, 
it is simple enough to obtain explicit results. 

In order to reduce ourselves to a base-only problem, we 
can follow the procedure outlined in the previous section. 
Attention must be payed to the fact that e indicates only 
the contribution of the electronic degrees of freedom to 
the full eigenvalue and the energy is given by e+ \k ^ x 2 . 
The Schrodinger equation for the base degrees of freedom 
can be obtained from the one in Ref. |ll| in two steps. 
First of all, the electron eigenvalue e must be substituted 
with /(e), defined in Eq. ([)]), a quantity that keeps track 
of the effects of the fibres on the electron dynamics. The 
second step amounts to re-normalize the wavefunction 
according to Eq. (||), i.e. - * ^[1 \J n {z) Hi - We then 
obtain, following Kabanov and Mashtakov |18| , 



(/(e) + 2AMeM| 2 )^ = *£^ 



+/ 



(12) 



It is now clear that the condition (JSj) introduces an ef- 
fective rescaling of the coupling constant A = g 2 /4kt. 
For n(e) < 1, the interaction term is less effective with 
respect to the base-only problem. Furthermore, every 
energy level can be mapped onto one of the base-only 
problem, with a coupling changing self-consistently with 
the energy value. 

As a consequence, the physical quantities involving ex- 
pectation values on a single eigenfunction will have simi- 
lar behaviors for the bundled structure and the base-only 
cases. On the other hand, we expect different properties 
for quantities which involve matrix elements between dif- 
ferent eigenfunctions, like e.g. spectral properties. 

In what follows, we choose the explicit case of the comb 
lattice, comparing our results with the ones of the one- 
dimensional case. The fibre is a semi-infinite linear chain 
and the expressions (||), (||) take the explicit form: 



on 



(e) = ^(e- S gn(eWe 2 -At 2 ) 



n(e) 



At 2 



(13) 
(14) 



Looking at eq. (|T|), we expect that, due to the effective 
reduction of the coupling constant, the crossover from a 
large to a small polaron, typical of the one dimensional 
case, will occur for a larger value of A. Moreover, the 
ground state properties of the one-dimensional case can 
be recovered from the ones of the comb lattice through a 
redefinition of the coupling and a rescaling of the wave- 
function: 



A -> An(e A ) 

ipi -* i>i/\Jn{e\) 

Xi — > Xi/ y/n(e\). 



(15) 



Here we emphasize the dependence on A of the eigenval- 
ues €\. 
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FIG. 1. Ground state energy for a 20 x 20 comb lattice (in 
units of t), as a function of the coupling constant A. 

We have also investigated if the rescaling of the eigen- 
value equation could induce a level crossing in the ground 
state at some value of the coupling, thus modifying the 
nature of the crossover with respect to the Id case. The 
first derivative of e(A) turns out to be proportional to the 
same quantity of the base-only case, through a continu- 
ous function of e, built up by Fqq(c) and its derivative. In 
this way we conclude that the properties of the ground 
state for a comb lattice, or generally for a bundled struc- 
ture whose fibres are semi-infinite linear chains, are the 
same of the ones of the base-only case. 

We performed numerical solution of the problem for 
finite, but quite large lattices, up to 20 x 20 sites. Hence- 
forth, energies will be expressed in units of the hopping 
amplitude t, and lengths in units of the lattice constant a. 
Our numerical data perfectly confirm the previous con- 
siderations. In figure [l] we plot the ground state energy 
as a function of the coupling constant A. It presents a 
continuous crossover between a large polaron state whose 
energy is only slightly modified with respect to the free 
value eo — — 2v2i to a strong coupling small polaron 
state, in which eo = —2ty/2X. This behavior is com- 
pletely analogous to the one-dimensional system |l8) . 

In order to estimate the crossover value A c , we plot 
the density-displacement correlation functions (\ipo\ Xi) 
for i = 0, 1, 2 (fig. ||), where i = indicates the site 
where the electron localizes Jlq ]. This quantity is a di- 
rect probe of the polaronic behavior, measuring the de- 
gree of correlation between the electron and the lattice 
deformations. All the non-local correlation functions de- 
cay rather sharply above a given value of A. Roughly at 
the same coupling value, the local correlation function 
abruptly changes the slope, signaling a crossover from a 
weak-coupling to a strong-coupling regime. This value 
can be associated to the crossover from a large polaron 
extended on several lattice site to a single-site small po- 
laron. 




X 

FIG. 2. Density-displacement correlation functions for a 
20 x 20 comb lattice for different distances of the sites from 
the density peak. The change of the slope of the local corre- 
lation function (bottom panel) is the signal of the polaronic 
crossover, and permits to extract the critical value A c = 0.925. 

We can therefore estimate the crossover coupling A c 
with the value where the first derivative of the local cor- 
relation function is maximum, or equivalently with the 
value for which the nearest-neighbours correlation func- 
tion starts to decrease. Using the first criterion we obtain 
an estimate A c = 0.925, which has to be compared with 
the one-dimensional value A c = 0.75. 

Notice that our comb lattice can be viewed as a first 
step beyond one-dimensionality. The crossover coupling 
for d = 2 is A ~ 1. One should anyway keep in mind that 
we are not considering el-ph interaction on the fibres. 

In Fig. [| we show how the Id quantities can be re- 
covered with the scaling relations (|lj): Once the local 
density-displacement correlation function (|^>o| &o) an d 
A are scaled respectively by 71 3 / 2 (£a) and n(c\), the curve 
lies on top of the one for the linear chain. 

As expected no qualitative difference with respect to 
the one-dimensional case has been found for the ground- 
state properties. 
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FIG. 3. The local density displacement correlation function 
for the linear chain is derived from the data by the comb 
lattice case with the scaling explained in the text (See Eq. 



A. Transport properties 

The study of the transport properties in the polaronic 
problem has been widely studied |2l],[l5|] , since it directly 
characterizes the onset of localization. 

One of the most striking features of bundled struc- 
tures is an anomalous diffusion for a classical walker mov- 
ing along the base ||. We expect quantum counterpart 
for this property, namely some peculiar behavior for the 
conductivity. Note that this expectation is not in con- 
tradiction with the general properties of the eigenstates 
described in the previous section, since the optical re- 
sponse is expressed in terms of matrix elements between 
different eigenstates. In fact the mapping (|l5| ) involves 
different normalizations, and therefore different effective 
coupling constants, for every energy level. 

The optical conductivity is given by: 



<t(uj) = D8(u) 



E 

n#0 



(16) 



where the coefficient of the zero-frequency delta function 
contribution is the so called Drude weight, which explic- 
itly characterizes the transport properties: a vanishing 
D is the signature of an insulating state ^2|, whereas a 
metal has a finite value of such a quantity. The Drude 
peak can be evaluated by means of the so called /- 
sumrule |23|: 



can be expected. We compute the Drude weight given by 
(|l6|), the total weight of the optical excitations, related 
to the kinetic energy by Eq. dl7|), and the integral over 
finite (non-zero) frequencies. Both the kinetic energy and 
the finite frequency optical conductivity show an anoma- 
lous behavior with respect to the one dimensional case as 
it is transparent from a comparison between the data in 
Fig. ^ and the ones for the one-dimensional case reported 
in the inset. 

For small values of the coupling constant the absolute 
value of the kinetic energy presents a small but visible 
enhancement with respect to its free value. This is con- 
firmed by a perturbative calculation for A < 1. In this 
limit, we can introduce a continuous approximation of 
eq. (H 



(a + 2t + 2\t\i}(x)\ 2 ) = -a 2 t 



2 d 2 ip(x) 



d x 2 

which has an exact solution of the form: 



ip(x) 



«( £ a) 



2a cosh(Xn(e\)x/a) 



(18) 



(19) 



In the one-dimensional case n(e\) — 1 and a just rep- 
resents the electronic contribution to the ground state 
energy, while for the comb lattice they are given by 
Eqs. (p"3|), and (|l4|). Thus the averaged kinetic en- 
ergy is Ekin = — 2t(l — X 2 t/6a 2 ) for the Id case, while 
E KIN ~ -V2t(l + X 2 t/8a 2 - A 4 /128a 4 ) on the comb lat- 
tice. This result can be easily understood by inspection 
of the particular geometry of the lattice. In the free case, 
the ground state wavefunction has a constant value on 
the backbone sites with exponentially decreasing tails on 
the teeth. Switching on the interaction, the backbone 
sites become energetically favorable and the particle is 
therefore recalled on the backbone, increasing the charge 
density on the substructure, and consequently of the 
kinetic energy. Further increasing A the localization ten- 
dency which leads to small polaron formation becomes 
effective. In this regime the usual behavior is recovered, 
the kinetic energy decreases its absolute value, while the 
e-ph interaction energy is substantially decreased. 

For small values of A the finite frequency contribution 
remains zero up to a given value of A, and continuously 
acquires a finite value at some value of A for both the 
structures. This can be attributed to finite size effects: 
a finite number of sites cuts off the available states with 
proper symmetry in (|l6|). 



a(uj)duj 



ire 



(H, 



kill I 



(17) 



and involves the calculation of the kinetic energy and of 
the integral over all finite frequencies of the optical con- 
ductivity. Quite naturally we only consider the conduc- 
tivity along the backbone, in which anomalous behavior 
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FIG. 4. Contributions to the DC conductivity along the 
backbone, compared with the Id case. In both cases we show 
the Drude weight D, the total sumrule for the optical con- 
ductivity (^), and the finite frequency optical weight. 

A striking feature emerges from Fig. ||. The coupling 
value for which the Drude weight vanishes (A ~ 0.72) 
is sizeably smaller than the crossover coupling extracted 
from the correlation functions (A c ~ 0.925), which in 
turn coincides with the kink in the kinetic energy curve. 
This is never the case in usual one-dimensional and two- 
dimensional lattices, where the loss of coherence and the 
small polaron formation essentially coincide, as can be 
seen in the inset, and as it is valid also for finite phonon 
frequencies. 

However, our results are not in contradiction with gen- 
eral criteria for small polaron formation fl5|l but they con- 
firm a crucial physical point. In fact the condition A > 1 
implies that the electronic energy to be overcome by the 
el-ph interaction energy for small polaron formation is 
not the coherent one but the total kinetic energy. These 
quantities coincide for a single particle on a conventional 
lattice, but they may differ for interacting systems or 
unordered structures. 

For A < A < A c the system is in a state in which 
the coherent transport is strongly suppressed while the 
kinetic energy is still large, substantially equal to its non 
interacting value. We attribute this contribution to local 
incoherent electron hopping between the teeth sites. 

This system represents, to our knowledge, the first one 
in which coherent motion is strongly inhibited, but there 
is a sizeable range of parameters in which no polaronic 
effects are found. This phenomenon represents the quan- 
tum heir of anomalous diffusion. It should be noted that, 
as the el-ph coupling is extended to the whole lattice this 



effect is expected to be emphasized. 

IV. CONCLUSIONS AND PERSPECTIVES 

The static and dynamical properties of an electron in- 
teracting with local classical oscillators on a comb lattice 
have been extensively analysed. Even though our work 
represents only the first attempt to tackle the problem 
of interacting fermions on unconventional structures, in- 
teresting and peculiar features have been discovered. On 
really general grounds, we have shown that the static 
properties (averages on a single state) of a general bun- 
dled structure can be easily obtained from the ones of 
the base-only problem by means of a rescaling of the pa- 
rameters. As an example, we have studied the small po- 
laron crossover on a comb lattice, comparing it with the 
one-dimensional case, explicitly confirming this property. 
Nevertheless, dynamical properties cannot be so simply 
recovered. A numerical study of the optical response on 
the comb lattice has highlighted two anomalous proper- 
ties. Contrarily to the one-dimensional case, for small 
e-ph couplings, the DC conductivity (measured by the 
Drude weight) is increased with respect to the bare non- 
interacting value. By increasing the coupling a rather 
surprising feature has been found, namely the coherent 
motion is strongly suppressed while a sizeable incoher- 
ent kinetic energy is still present. The loss of coherence 
is therefore not due to the self-trapping associated with 
the polaron crossover (that occurs for larger e-ph cou- 
pling, when also the incoherent kinetic energy is sup- 
pressed), but to the drastic effects of a unconventional 
geometry on the electronic properties. The occurrence of 
such an effect in classical systems has been demonstrated 
and widely discussed within the framework of the ran- 
dom walk. We note that our approach is not completely 
the analogous of the random walk problem: the absence 
of e-ph interaction on the teeth sites implies that, in the 
classical analogous, the motion of the walker on these 
sites could be deterministic. The quantum counterpart 
of anomalous diffusion driven by the geometry would be 
achieved only if the e-ph interaction was extended to the 
whole lattice. 
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